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Abstract. We study the dynamics of the thermal momentum distribution 
function for an interacting, homogenous Fermi gas on 7? in the presence of 
an external weak static random potential, where the pair interactions between 
the fcrmions arc modeled in dynamical Hartrcc-Fock theory. We determine the 
Boltzmann limits associated to different scaling regimes defined by the size of 
the random potential, and the strength of the fermion interactions. 



1. Introduction 

We study the dynamics of an interacting, homogenous Fermi gas on Z 3 in a 
static, weakly disordered random medium, where the pair interactions between the 
fcrmions are modeled in dynamical Hartree-Fock theory. An observable of consid- 
erable importance is the momentum distribution function at positive temperature, 
and we are interested in its dynamics for time scales that are associated to kinetic 
scaling limits of Boltzmann type. A main motivation is to understand the trend to 
equilibrium in such systems, and to control the interplay between the influence of 
the static randomness, and nonlinear self-interactions of the particles. We derive 
Boltzmann limits for the thermal momentum distribution function, depending on 
different scaling ratios between the random potential, and the strength of the pair 
interactions between the fermions. 

The model in discussion describes a gas of fermions in a finite box A^ := 
[— ^, j] 3 n T? of side length L » 1 with periodic boundary conditions, and as- 
sociated dual lattice A^ := Al/L C T 3 ; we will eventually take the thermo- 
dynamic limit L — > oo. On the fermionic Fock space of scalar electrons, $ = 
C © ©„>i Ai ^ 2 (Al)j we introduce creation- and annihilation operators a+, a y , 
for x, y € A^, satisfying the canonical anticommutation relations {a~£,a y } :— 
a x a y + a y a t = L 3 5 Pi q where S p . q is the Kronecker delta, and {a|,a$,} = for 
a' = a or a+. We denote the Fock vacuum by fi = (1, 0, 0, . . . ) e it is annihilated 
by all annihilation operators, a x il = for all x C A^. 

Letting 21 denote the C*-algebra of bounded operators on we consider a time- 
dependent state p t on 21 determined by 

id tPt (A) = Pt ([H(t), A]) (1.1) 

for A g 21, with a translation invariant initial condition p . We assume that p is 
number conserving p ([A,N}) = for all A £ 21, where N := J2 x eA L a t a x denotes 
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the particle number operator. We study the dynamics of p t determined by the 
time-dependent Hamiltonian 

H(t) = H o + nV u + XW(t) , (1.2) 
where Hq := / dpE(p) a+a p is the second quantization of the centered nearest 
neighbor Laplacian, with E(p) = 2 X^j=i cos ( ^Pj ) denoting its symbol. The 
interaction of the fermion gas with the static random background potential is de- 
scribed by the operator V u := J2 x eA L u ^ a t a x where {w x } lg A L is a field of i.i.d. 
centered, normalized, Gaussain random variables. The small parameter < r\ <C 1 
determines the strength of the disorder. 

The operator W(t) accounts for fermion pair interactions in dynamical Hartree- 
Fock theory, 

W(t) := v(x-y)[E[p t (a+a x )]a+a y - E[p t (a+a x )]ata y ] . (1.3) 

Here, v denotes a pair potential, where ||u|| i3 -3/2+<x( T 3) < C is assumed for a > 
arbitrary but fixed. Notably, the unknown quantity E[ p t (a x a x ) ] itself appears in 
W(t). 

Our main interest is to study the dynamics of E[p t (a+a x )], the average of 
the pair correlation function, which is determined by the self-consistent nonlinear 
evolution equation 

id t E[p t (ata y )} = E[pt([H{t) , a+ay])] (1.4) 

with initial condition p (a x r a y ). In particular, we derive Boltzmann equations in 
kinetic scaling limits of the above model, for scaling regimes defined by different 
ratios between rj and A. 



The relevant scaling relations in the system can be understood with the help 
of the following heuristics. The assumption of vanishing mean implies that the 
average effect of the random potential on the dynamics of p t in a time interval [0, t] 
is proportional to its variance, of size 0(r/ 2 t). This suggests that the strength of 
the pair interactions between fermions, and the interactions of each fermion with 
the random potential, are comparable if A = 0(ij 2 ). Accordingly, we distinguish 
the following scaling regimes, for which we derive the associated Boltzmann limits: 

• The scaling regime A = 0(ry 2 ), where the interactions between pairs of 
fermions, and of each fermion with the random potential are comparable. 
For any T > 0, for test functions /, g, and T/t — rj 2 , we prove 



hm km E[ PT/J)2 (a+ (/)«(<?))] - / dp f(p) g(p) F T (p) , (1.5) 

rj— >0 L— s-oo J 

where a + (/) = ^2 x f(x)a x = (a(/))* (adjoint), and Ft(p) satisfies the 
linear Boltzmann equation 

d T F T (p) = 2n J duS(E(u) - E(p) )(F T (u) - F T (p) ) (1.6) 

with initial condition Fo = po e H^ +a (T 3 ) for some a > 0. The fact 
that the resulting Boltzmann equation is linear follows from complicated 
phase cancellations due to translation invariance. While the microscopic 
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dynamics is nonlinear, we prove that its kinetic scaling limit, as n — > 0, is 
described by a linear Boltzmann equation. The scattering kernel accounts 
for elastic collisions preserving the kinetic energy. This result remains valid 
in the regime A = o(n 2 ). In the case A = 0, it reduces to the one proven in 

D3B- 

In the regime rf = o(A) and (T, X) = (At, Ax), we prove that the kinetic 
scaling limit A — > is stationary. 

In the regime where A > is independent of rj, and for the scaling (T, X) = 
(r] 2 t, r] 2 x), we characterize the stationary states; those are given by solutions 



of an implicit equation of the form (3.4) with zero on the l.h.s., but where 



the delta distribution enforces conservation of a renormalized energy per 
particle. Accordingly, the stationary states are supported on level surfaces 
of a renormalized kinetic energy function, determined by a nonlinear fixed 
point equation. A derivation of non-stationary solutions in this kinetic 
scaling limit is an interesting open problem. 

Our work significantly extends (TUJ which addresses the Boltzmann limit for a 
homogenous free Fermi gas in a random medium. The proofs given in |10| employ 
techniques developed [3 [HI H2 developed for the derivation of Boltzmann equa- 
tions from the quantum dynamics of a single electron in a weak random potential; 
see also [22, 26 . In the landmark works [T3 1 fl4| . [To] . this analysis has been extended 
to diffusive time scales. We also refer to [TJ [HI [TTJ US for related works. 

For the proof of our results, we represent the average momentum density E[/i t ] 
in integral form, as an expansion organized in terms of Feynman diagrams. Our 
overall strategy parallels the approach in [T7J [HI EH EU 03 1211 HH] , but the techniques 
developed in those works (for linear models) do not carry over directly because of 
the nonlinear self-interaction of the fermion field. In particular, resolvent meth- 
ods which underlie the analysis those works are not available here. Instead, our 
approach strongly uses stationary phase techniques, in order to control the combi- 
nation of Feynman graph expansion techniques with such nonlinearities. 

For the related problem of the derivation of dynamical Hartree-Fock equations 
from a fermion gas, see for instance [3]. We note that the derivation of macro- 
scopic transport equations from the quantum dynamics of Fermi gases without any 
simplifying assumptions on the intcrparticle interactions (and without random po- 
tential) is a prominent and very challenging open problem in this research field; see 
for instance [SJ [TCI HH H31 [27] . For some very interesting recent progress relevant 
related to this issue, see |2"4"1. 
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2. Definition of the model 

We consider a gas of fermions in a finite box A^, := [— ^, D 1? of side length 
L » 1, with periodic boundary conditions. We denote the associated dual lattice 
by A^ := Al/L C T 3 . We assume that L is much larger than any other significant 
length scale of the system, which will depend upon the case under consideration. 
We will eventually take the thermodynamic limit L — > oo. 

We denote the Fourier transform by 

f(p) := £ e-^/Or), (2.1) 

x£A L 

where p e A^ , and the inverse transform by 

g y (x) = J dpe 2 ^ x g(p). (2.2) 
For brevity, we are using the notation 

dpf(p) = ^ /(P)' ( 2 - 3 ) 



which recovers its usual meaning in the thermodynamic limit L — > oo. 
We will use the notation 

S(k) := LH k , (2.4) 

where 

f 1 if fc = 

dfe ~ \ otherwise (Z ' b) 
denotes the Kronecker delta on the momentum lattice A^ (mod T 3 ). 
We denote the fermionic Fock space of scalar electrons by 

3 = 03™, (2-6) 

n>0 

where 

n 

do = C , $ n = f\e 2 (A L ) , n> 1. (2.7) 
i 

We introduce creation- and annihilation operators a+, a q , for p, q e A* L , satisfying 
the canonical anticommutation relations 

a+ a q + a q a+ = *(p otherwise. (2 " 8) 

There is a unique unit ray 1! = (1,0,0,...) £ J referred to as the Fock vacuum, 
which is annihilated by all annihilation operators, a p f2 = for all p G A* L . 

Let 21 denote the C* -algebra of bounded operators on We consider a time- 
dependent state p t on 21 determined by 

id t p t (A) = p t ([H(t), A}) (2.9) 
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for A € 21, with a translation invariant initial condition pq. We assume that po is 
number conserving; that is, 

p ([A,N})=0 (2.10) 

for all A e 21, where 

JV := £ a+a x (2.11) 

denotes the particle number operator. 

The dynamics of p t shall be determined by the time-dependent Hamiltonian 

H(t) = H o + + \W(t) , (2.12) 
where the right hand side is defined as follows. The operator 

H := J d P E(p)a+a p (2.13) 

is the second quantization of the centered nearest neighbor Laplacian (Af)(x) = 
J2\ y -x\=i f(v) 011 The symbol of A is given by 

3 

E(p) = 2^cos(27r Pj ), (2.14) 

corresponding to the kinetic energy of a single electron. The interaction of the 
fermion gas with the static random background potential is described by the oper- 
ator 

V u := a * ■ ( 2 - 15 ) 

xeA L 

We assume {co x } xe \ L to be a field of i.i.d. random variables which is centered, 
normalized, and Gaussian, 

EK] = , E[oj 2 x ] = 1, (2.16) 

for x <G Al. The small parameter 0<f)<l controls the strength of the disorder. 

The operator W{t) accounts for fermion pair interactions in dynamical Hartree- 
Fock theory, 

W(t) := v(x-y)[®[pt(a+a x )]a+a y - E[p t ( a+a, ) ] a+a y ] . (2.17) 

x,yeA L 

Here, v denotes a A^-periodic pair potential, for which we assume that 

|| v\\ H3/ 2+, {T3) < C (2.18) 

for a > arbitrary but fixed. Notably, we make no assumption on the sign of v. 
We note that the unknown quantity E[p t (a+a x ) ] appears in W(t). 

We are interested in the dynamics of the average of the pair correlation function, 

E[p t (a+a x )], (2.19) 
which is determined by the self-consistent nonlinear evolution equation 

id t E[pt(ata y )} - E[ Pt ([H(t),4a y })} (2.20) 
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with initial condition po(a x a y ). Its Fourier transform is diagonal in momentum 
space, 

E[p t (a+a q )} = 5(p-q)^E[p t (a+a p )], (2.21) 

because E{p t (a,ya x )] is translation invariant. This follows from the homogeneity 
of the randomness. 

We remark that for fermions, 

< ^p (a+a p ) < 1, (2.22) 

since ||ap + ^|| = L d l 2 in operator norm, Vp € A^. The expected occupation density 
of the momentum p in the lattice is given by 

VH(p) := pE^o+Op)]. (2.23) 

The dynamical Hartree-Fock interaction can be written as 

W(t) = ^p (N)(Y,v(x))N - J dp(v * p t )(jp)a+a p 

=: W dir (t) + W ex (t) (2.24) 

where, following standard terminology, Wdi r {t) denotes the direct, and W ex (t) the 
exchange term. It is clear that H^it) is particle number conserving, 

[H u (t),N] = , VieK. (2.25) 

Since po is translation invariant and number conserving, we conclude that whenever 
[A, N] = holds for an operator A, it follows that 

id t p t (A) = p t ([H ex (t),A}) (2.26) 

where the operator 

H ex (t) := H + ryVL + XW ex (t) (2.27) 
contains only the exchange term of W(t). 
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3. Statement of the main results 



In order to determine the dynamics of the average momentum distribution func- 
tion p t defined in (2.23), we consider 

dpWj g( P ) Mp) = E[ Pt (a+(f)a(g))] 

= E[p {U* t a + {f)a{g)U t )} (3.1) 

for a translation invariant and particle number conserving initial state po, where / 
and g are test functions. The linear operator lAt denotes the unitary flow generated 
by H ex (t). It satisfies Uq — 1, and notably depends on p si s 6 [0,4]. 



Accordingly, we make the key observation that (3.1 ) is a fixed point equation for 
Pt, tested against /, g. The right hand side of (3.1) is a complicated nonlinear 



functional of p s which we will discuss in detail in Section 4.1 



We introduce macroscopic variables (T, X), related to the microscopic variables 
(t,x) by 



(T,X) = (Ct.Cz): 



(3.2) 



with ( > a real parameter. We will study kinetic scaling limits associated to 
different scaling ratios between £, n and A. 

As stated in the introduction, the random potential has an average effect on 
the dynamics of p t by an amount proportional to its variance, 0(rj 2 t), in the time 
interval [0,4]. Since the strength of the fermion pair interactions is O(X), both 
effects are comparable if A = 0(r] 2 ). This implies that the relevant scaling regimes 
of the system are determined by those addressed below, in Theorems 3.1 3.2 and 
Theorem 13.41 



Theorem 3.1. Assume that A < 0(r] 2 ). Then, for any fixed, finite T > 0, and 
any choice of test functions f , g, 



(3.3) 



Urn lim E[p T/v ,(a+(f)a(g))) = / dp f(p) g(p) F T (p) 
holds, where Ft (p) satisfies the linear Boltzmann equation 

d T F T {p) = 2n J duS(E(u) - E(p)) (F T (u) - F T (p)) 
with initial condition Fq = po 6 iJ2 +<T (T 3 ) for some a > 0. 



(3.4) 



We note that the linear Boltzmann equation (3.4) can be explicitly solved. The 
solution is given by 



F T (P) = Foo(p) + e~ Tm ^ 
where 



2tt 



du5(E(u) - E(p) ) ( F (u) - F (p) ) , (3.5) 



m(p) 

m(p) := 2tt / duS(E(u) - E{p) ) 



(3.6) 
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and 

2-7T f 

:= —r -r / du$(£(u)-.E(p))Fo(tt). (3.7) 
m(p) 7 

As an important example, we note that the following is obtained if the initial 
state po is given by the Gibbs state for a non-interacting fermion gas (with inverse 
temperature j3 and chemical potential p), 

Po (A) = -^Tr(e-^-^U), (3.8) 

where ^ := Tr( e-P( T ~^ N ) ). The associated momentum distribution function is 
given by the Fermi-Dirac distribution 

2™o h Po( <%> ) = 1 + Je {p) -») > ( 3 - 9 ) 



which is a stationary solution of the linear Boltzmann equation (3.4 1, for all (3 > 0. 
This result remains true in the zero temperature limit (3 — > oo where, in the weak 
sense, 

X[E( P ) < fi] (3.10) 



1 + e P(E(p)-n) 

(see also [TO]). 

In the case ij 2 — o\(l) and T = Xt, we find the following kinetic scaling limit. 

Theorem 3.2. Assume that rf — 0(X 1+S ) for S > arbitrary. Then, for any 
fixed, finite T > 0, and any choice of test functions f , g, 

lwxlmiE[p T/x (a+(f)a(g))}= fdpf(pjg(p)F T (p), (3.11) 

A— fO L— >oo J 

and 

d T F T (p)=0, (3.12) 
for Fo = po 6 i/5 +cr (T 3 ) /or some <r > Accordingly, Ft = Fo is stationary. 

Finally, we prove a partial result that highlights some interesting aspects about 
the problem of determining the kinetic scaling limit determined by T = rj 2 t and 
77 — 5- 0, with A small but independent of 77. That is, we are considering, for A = O(l), 
the rescaled, formal fixed point equation 



/ 



dpf(p)g(p)PT/^(p) = G (L} [p,{»);vA; T ;f,g} 

:= E[p T/ri 2(a+(f)a(g))} (3.13) 

for p,, (•). The existence and uniqueness of solutions for this fixed point equation 
is currently an open problem. Below, we will make the assumption that there exist 
limiting stationary solutions, and determine a their form under this hypothesis. 

We base our discussion on the following hypotheses for the case A = O(l): 
(HI) There exist solutions F^(T) :— lim^^co pTjrj 2 of (3.13), such that the 



limit w - lim^o F^(T) =: F(T) = F(0) exists and is stationary. 
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(H2) The stationary fixed point solution in (HI) satisfies 

F{T) = lim lim G^[F ( ^ ; m X;T; f, g] 

77— J>0 L — >oo 

= lim l lim g^[F; m \;T;f,g}. (3.14) 

rj— fO L— >oo 

The first equality sign here is equivalent to (HI), while the second equality 
sign accounts for the assumption that can be replaced by the limiting 
fixed point F before letting r\ — > 0, to produce the same result. 

We remark that based on the analysis given in this paper, we are able to prove 
hypothesis (H2) if F^) = F + 0{rf ). Error bounds of order 0(r] 2 ) require more 
precise estimates of "crossing" and "nesting" terms in the Feynman graph expansion 
than considered in this paper, but are available from [T31 Q31 H51 US] ■ We will not 
further pursue this issue in the work at hand. 

Proposition 3.3. Let A be small but independent of r\, and assume that F £ 
L°°(T 3 ) independent oft. Then, the thermodynamic limit 

g[F; m X;T;f,g] := lim G^[F; V ; A; T; f,g] (3.15) 

L— >-oo 

exists. 

The proof of this proposition follows directly from results established in HI 
[TOl fm , and will not be reiterated here. 



Theorem 3.4. Assume that A < O v (l), and let 

E x (u) := E{u) - \(v*F){u). (3.16) 
We assume that F G _L°°(T 3 ) admits the bounds 

sup { dp 1 -, sup I da 1 -<Clogi (3.17) 

a J \Fj\{q) — a — ie\ g J \h\(q) — a — ie\ e 



nil 



d 



sup sup / dq dp 



«f «£pJ ~ \E\(q) - ax - ie\ \E\(p) - a 2 - ie\ \E\(p ± q + u) - a 3 - ie\ 

< e- b (3.18) 

for some < b < 1. 

Then, F satisfies 

[ dpJ(p)g{p)F{p) = lim G[F;ri;X;T;f,g], (3.19) 

independent of T, if and only if it satisfies 

F( P ) = Mob) = / du6(E x (u)-E x (p))F(u), (3.20) 
m\{p) J 

where 

fh x {p) := 2tt J du8{E x (u) - E x {p) ) (3.21) 
is the (normalized) measure of the level surface of Ex for the value E x (p). 
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Remark 3.5. The following comments refer to Theorem\3.4 



(1) The solution of (3.16) corresponds to a renormalized kinetic energy which 
is shifted by the average interaction energy for fermion pairs. 



(3) 



(2) The fixed point equation (3.20) for F shows that the stationary kinetic limits 



of /it are concentrated and equidistributed on level surfaces of the renormal- 
ized kinetic energy function E\{ ■ ). 



The bounds (3.11) and \ 3. Mfy correspond to the "crossing estimates" in 
ensure sufficient non- degeneracy of the renormalized 
energy level surfaces so that the Feynman graph expansions introduced below 
are convergent. However, they do not seem sufficient to prove hypothesis 
(H2) under the assumption that (HI) holds. 



(4) We note that if A < o^(l), the stationary solutions found in Theorem 3.4 



reduce to those of the linear Boltzmann equation derived in Theorem \3.1 
see (3~l\). 
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4. Feynman graphs and amplitudes 

In this section, we set up the Feynman graph expansions underlying our proofs 
of Theorems [O [3~2l and Theorem IOI 



4.1. Duhamel expansion. Let Ut denote the unitary flow generated by H ex (t), 
determined by 

id t U t = H ex (t)U t and U = 1 . (4.1) 

It then follows that 

Pt(A) = po(U;AU t ), (4.2) 

and that 

pt(a+(f) a(g) ) = p (a+{f, t) a(g, t) ) . (4.3) 
The Heisenberg evolution of the creation- and annihilation operators is determined 

by 

a{f,t) :=W t a{j)U t , (4.4) 

with 

a(f) = J dpf(p)a p , a+(f,t) = (o(/,t))*. (4.5) 

It suffices to discuss the annihilation operators a(f, t). Because H ex {t) is bilinear in 
a + and a, it follows that a(f, t) is a linear superposition of annihilation operators. 
Therefore, there exists a function f t such that 

a(f,t)=a(f t ), (4.6) 

satisfying 

id t a(f t ) - [H ex (t),a(f t )} 

dpft{p)E{p)a p + T] J dp J duf t (p)Q(u- p)a u 

-A J dp(v* nt)(p)ft(p)a p , (4.7) 

with initial condition 

o(/,0) = a(/ ) = a(/). (4.8) 
We conclude that f t is the solution of the 1-particle random Schrodinger equation 
id t ftip) = E(p)f t (p) + v (Q*f t )(p) - \ (v* iH)(p)f f (p) (4.9) 
with initial condition 

fo = f. (4.10) 
Here, Q(u) — J2 X e 2vlux uj x , and v is the fermion pair interaction potential. 

Noting that the Hamiltonian H ex (t) itself depends on the unknown quantity /it, 



we determine \i t by writing the fixed point equation (3.1| in integral form, as an 
expansion in powers of rj. 
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For arbitrary test functions / and <?, we consider the pair correlation function 
Pt(a + (f)a(g) ) = Po(a+(f t ) a(g t ) ) 



dp dq po ( a+ a q ) f t (p) g t {q) 



dpJ{p)ft{p) gt(p) 



(4.11) 



where the state p t equals the one in the definition of /i t , (2.23). Passing to the last 

(4.12) 



line, we have used the momentum conservation condition 

Po(«p a q ) = J{p)5{p-q) 
obtained from the translation invariance of the initial state po, where 



< J(p) = -^p (a+a p ) < 1, 



(4.13) 



as noted before, in (2.22 1. 



The solution f t of (4.9), (4.10), satisfies the Duhamel (respectively, variation of 
constants) formula 



ft(p) = U , t (p)f(p) + ir) dsU Stt (p)(u)*f s )(p) 



with 



(4.14) 
(4.15) 



U a>t (p) := ^^f^l-^W), 
where we treat 

K s (u) := (v*p s )(u) (4.16) 

as an external (a priori bounded) source term. We note that Uo^{p)f{p) solves 
(4.9) for ij — (no random potential) with initial condition (4.10). 



Let N £ N, which remains to be optimized. The AT- fold iterate of (4.14) produces 
the truncated Duhamel expansion with remainder term, 



f _ A<N) f (>N) 

h — it + ft 



where 



N 



f(<N) 



(«) 



and ft >N ^ is the Duhamel remainder term of order N. We define 



so 



(4.17) 



(4.18) 



(4.19) 



t-i := , tj 
for j — 0, . . . , n, and 

r n t ■ 

n(k ,...,k n ;z) := / ^so-.-ds^Re-^^e^V-^*'^'^)^) 



for z € 
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The n-th order term in the Duhamel expansion is given by 



jf := (ir,r I dt 







/ dt n -- 




Jo 





dti / dk ~- dk n 8{p - k ) (4.21) 



It 74 

[ n u **-i* (%) ] [ n ^ - ^-i) ] /(*») • 

3=0 1=1 

Expressed in terms of the time increments Sj :— tj — tj-±, 

ti n \p) = (^7 7 )" I ds ■ ■ ■ ds n 5(t ~ s j) / dk ° ■ ■ ■ 6 <P ~ fc ") ( 4 ' 22 ) 
J j=o J 

n t n 

Expressing the delta distribution 5(t — Y^j=o s j) m t erms of its Fourier transform, 
we find 

ft n) (P) = (iv) n e £ * / da e~ Ua J dk ■ ■ ■ dk n S(p - k ) 

n 

JJcoikj -kj-i) f(k n ). (4.23) 

3=1 

The above three equivalent expressions for (p) have different advantages which 
we will make use of. 



The Duhamel remainder term of order N is given by 



/<>"> = «,/ dsU a MVfW 



We choose 



t 



(4.24) 
(4.25) 



so that the factor e in (4.23) remains bounded for all t. 



Substituting the truncated Duhamel expansion for a + (ft), a(g t ) in (4.111, one 
obtains 



N+l 



Pi 



(a+(f)a(g)) = p (a+(f t ) a(g t )) = £ p?> % \f,g) 



n,n— 



where 



if n, n < N, and 



pl n ' N+1) (f,9) := Po( a + (/ t (n) )a(ff t (>JV) )), 

(N+l.n), f , ._ / + / f{>N)\ i (n) s 



Pi ■ ' US) :=Po(a + (/r V >(ffr)) 
if n, n < AT. Moreover, 

^ +W1> (/,S) := Po(« + (/ t (>JV V(<? t (>Ar) ))- 



(4.26) 

(4.27) 

(4.28) 
(4.29) 

(4.30) 
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In particular, all Duhamel terms indexed by n, n < N depend on Q like polynomials. 
Accordingly, 



7rer„ 



t rt 2 rt _ 

dt„ ■ ■ ■ / dti / dig 



dtx 



du --- du 2fl +i f(uo) g(u 2 n+i) J (un) S(u n ~ u n +i) 

n 2n+l 



3=0 



j=n+l 



(4.31) 



2n+l 



E 



[ n ~ u i-^ n ~ 

j=n+2 



3 = 1 



and using (4.23), this is equivalent to 



E[pi"' S) (/,g)] = V m e^ £ dad5e« a -*> 



du --- du 2fl +l f(uo) g(u 2 n+l) J (tin) *(Un - tt n+ i) 

7^(tt , . . . , ti„; a + te) 7£(ti n+ i, . . . , tt 2 „+i; a - ie) (4.32) 

n 2n+l 
E [ JJ D(<!Zj - U 3 -_l) - Uj-l) 

3=1 j=n+2 



where t_i,t_i := in (4.31) 



4.2. Graph expansion. By assumption, {ui x } is a centered, i.i.d., Gaussian ran- 
dom field. Accordingly, we may explicitly determine the correlations of the random 
potential in the expressions (4.31), ( 4.32[ ). The expectation of any product of even 
degree n € 2N is equal to the sum of all possible products of pair correlations of 
the same degree, 



E[jj "(«,■-«,■-!)] 



(4.33) 



3=1 



pairings (ti,£'i) i=l 

The sum extends over all possible pairings (4,^) S {1, . . . ,n} 2 with £j ^ z = 
1, . . . , ^, where every element of {1, . . . , n} appears in precisely one pairing. This 
expansion is often referred to as Wick's theorem. Expectations of products of uj of 
odd degree are identically zero. 

To organize the terms in these sums of products of pair correlations, we introduce 
Feynman graphs. The set of Feynman graphs T n ^, with n = n + n £ 2N, is given 
as follows; see also [51 [TU1 HZ] : 
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] — 6- 



We consider two horizontal solid lines, which we refer to as particle lines, 
joined by a distinguished vertex which we refer to as the po~ ver tex (corre- 
sponding to the term Po(oi„ a u n+1 )• See Figure 1 for an example. 
On the line on its left, we introduce n vertices, and on the line on its right, 
we insert n vertices. We refer to those vertices as interaction vertices, and 
enumerate them from 1 to 2fi starting from the left. 

The edges between the interaction vertices are referred to as propagator 
lines. We label them by the momentum variables uq, W2n+i, increasingly 
indexed starting from the left. To the j-th propagator line, we associate the 
propagator Ut j _ 1 ,t j (uj) if < j < n, and Uj._ l (uj) if n + 1 < j < 2n + 1 
(with reference to the expression ( |4.31 1). 



To the £-th interaction vertex (adjacent to the edges labeled by u^-i and 
ui), we associate the random potential u)(u£ — ue~±), where 1 < £ < 2n + 1. 

n i "---! , 2n 

— □ □ □ ■ tiz 6 □ □ □ □ □ 

T 1 n+1 j i ■ p 



Figure 1. A contraction graph. 



A contraction graph associated to the above pair of particle lines joined by the 
Po-vertex, and decorated by n + n interaction vertices, is the graph obtained by 
pairwise connecting interaction vertices by contraction lines. We denote the set of 
all such contraction graphs by r rlj j;; it contains 

\T nfi \ = (2ft - l)(2n - 3) • --3 • 1 = ^ = 0(n!) (4.34) 

n!2" 

elements. 

If in a given graph ir £ T n %, the £-th and the £'-th vertex are joined by a 
contraction line, we write 

£ ~„ £' , (4.35) 
and we associate the delta distribution 

5{ui - ue-i - (ui> - ui'-i)) = E[Q(u£ - ui-i )Q(ue> - uy-i )] (4.36) 
to this contraction line. 

We classify Feynman graphs as follows; see also [8l HT]: 

- A subgraph consisting of one propagator line adjacent to a pair of vertices 
£ and £ + 1, and a contraction line connecting them, i.e., £ ^ £ + 1, where 
both £, I + 1 are either < n or > n + 1, is called an immediate recollision. 

- The graph 7r € T n ^ n (i.e., n = n = n) with £ ^ 2n — £ for all £ = 1, . . . ,n, 
is called a basic ladder diagram. The contraction lines are called rungs of 
the ladder. We note that a rung contraction always has the form £ ^ n £' 
with £ < n and £' > n + 1. Moreover, in a basic ladder diagram one always 
has that if l\ ~ n £[ and £ 2 — n £' 2 with £ x < £ 2 , then f 2 < £[. 
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- A diagram tt £ r„^ is called a decorated ladder if any contraction is cither 
an immediate recollision, or a rung contraction ij ~ w £'j with tj < n and 
^ > n for j = 1, . . . , k, and £\ < ■ ■ ■ < £^, £'\ > ■ ■ ■ > £' k - Evidently, a basic 
ladder diagram is the special case of a decorated ladder which contains no 
immediate recollisions (so that necessarily, n = n). 

- A diagram tt £ T n _„ is called crossing if there is a pair of contractions 
£ ~tt j f, with £ < £' and j < f , such that £ < j. 

- A diagram tt £ T n „ is called nesting if there is a subdiagram with i ~„. 
£ + 2k, with k > 1, and either I > n + 1 or £ + 2k < n, with j ^ j + 1 for 
i j=£+l,£ + 3,...,£ + 2fc— 1. The latter corresponds to a progression of 

— 1 immediate recollisions. 

We note that any diagram that is not a decorated ladder contains at least a crossing 
or a nesting subdiagram. 



4.3. Feynman amplitudes. To every Feynman graph tt £ T n .n we associtate its 
Feynman amplitude, as follows. 



To start with, we recall from (4.201 



U{u , . . . ,u n ;a + ie) (4.37) 



where = sq + • • • + Sj for j > 0, and t_i =0. Moreover, we recall that 



K s (p) = (v*fi s )(p) (4.38) 

where 

< n s {p) < 1 (4-39) 
holds uniformly in s and p. As a consequence, 

II ''II,., , = / dpi Y^v{x)e- 2 ^ x 
Jt 3 1 „ 



k s||l~(t 3 ) S ||w||li(t 3 ) 



< Vol{T 3 } || ( • v || £2(z3) || ( • )- 3 / 2 - ||, 2(z3) 

< C || U 11^3/2 + =- 

< C" (4.40) 



uniformly in s G K + . Here, we have recalled the property (2.18) satisfied by the 



fermion pair interaction potential v, for constants C, C that depend on a > 0. 
Given tt £ T n ^, we define 

^({%}j=o X ) := II (4-41) 

where every contraction line in 7r corresponds to one of the factors on the rhs, with- 
out any repetitions. Clearly, S v determines the momentum conservation conditions 
on the graph n. 
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Definition 4.1. The Feynman amplitude associated to the graph n € r„,^ is defined 
by 

:= ,»e»/ M e*-> (4.42) 

du ■ ■ ■ du 2 n+i f(uo) g(u 2 n+i) J (it„) S(u n - U n+ i) 
M {uj^lo 1 ) 1l(uo, ...,u n ;a + ie) U(u n+ i, . . . , u 2 r+i;5 - ie) • 

Our choice of e will be e = j. 

Since {w x } are i.i.d. centered Gaussian, 

E[Y[u(u t -ut-i)] = im^iui-Ui^uiiuj-Uj^)] (4.43) 

equals the sum of all possible products of pair correlations 

E[£D( u)lj(u')] = S(u + u') (4.44) 
(Wick's theorem). Accordingly, 

E[^"' S) (/,S)] = E Amp^/,^) (4.45) 

is the sum of Feynman amplitudes of all Feynman graphs n € r„ f;. 
As a consequence of translation invariance of po, we have that 

Pa( a t n a u n+ i ) = ^Mn - «n+l ) , (4.46) 

as we recall from (4.12). Moreover, translation invariance also implies overall mo- 
mentum conservation, that is, 

u - u 2n +i = , (4.47) 
which one easily verifies by summing up the arguments of all delta distributions. 

Accordingly, we arrive at the expansion 

N+l 

®[Pt(f,9)] = Amp^(/, 5 ;e;7 ? ) 

n,n— 7TEr n) $5 
N 

+ Y,(np^ N+1) (f,9)] + np^ +1 ' n \f,9)}) 

n=0 

+ n P ^ +1 ' N+1 \f,9)] (4.48) 
where the hrst term on the rhs is entirely expressed in terms of Feynman graphs 
and Feynman amplitudes. The terms on the second and third line on the r.h.s. 
involve the Duhamel remainder term, and will be shown only to contribute to a 
small error. 



18 



T. CHEN AND I. RODNIANSKI 



5. Proof of Theorem I3JJ I. Bounds on error terms 



In order to prove the Boltzmann limit stated in Theorem |3.1| we separate the 



main terms in the expression (4.48) from the error terms. We subsequently show 



that the main terms converge to a solution of the Boltzmann equation, while the 
error terms tend to zero. The Feynman graphs associated to the main term cor- 
respond to those appearing in the works [T71 O [10]. However, many aspects of 
the approach developed in those works (for the weakly disordered Anderson model, 
which is linear) are not suitable for the problem at hand. The main issue is the 
presence of the phase A J^ +s K s i(u)ds' in (4.37), which depends on the unknown 



quantity /it itself. In this section, we introduce a main tool, given in Lemmata 5.1 



and|5.2[ that enables us to control the nonlinear self-interaction of the fcrmion field. 



In a first step, we prove an estimate that will serve as a substitute for resolvent 
estimates. The latter were abundantly used in |17[ 15] . but due to the nonlinear 
self-interactions of the fcrmion field, they are not available here. 

Lemma 5.1. Let e = j <C 1. Then, there exists a constant C < oo independent of 
n, X, e such that 

A \ "+ 1 A 1 

(i \E(uj) - a\ + e 



\n(u ,...,u n ;a + ie)\ < C n+1 ( 1 + - J J] 
for all n G N. 



(5.1) 



3=0 



Lemma 5.1 follows from Lemma 5.2 below. We note that in our proof below, 
we allow for at most one integration by parts with respect to the time variable 
s, which requires no smoothness assumption k s (u). Iterating integration by parts 
with respect to s would easily produce the asserted result, but under the assumption 
that k s is C™-smooth in s. However, we will only rely on the a priori boundedness 
of k s (u), which is a consequence of the Fermi statistics satisfied by the quantum 
field, and will not assume any smoothness with respect to s. 



Lemma 5.2. Assume (4-4ty - Then, uniformly in r > 0, 

/ (u)ds' 



dse- ls{E{u) - a - ie) e~ iX K 



7 + 



< 



A . 
1 + - 
e 



C 



\E(u) - a\ + e 



where E{u) is the symbol of the nearest neighbor Laplacian on 1? . 



(5.2) 



Proof. We define 



Kt,t+s(u) 



t+s 



ds' K s ' (u) . 



Clearly, |K t]t+s (u)| < Co, uniformly in t and s > 0. 



The integral on the left hand side of (5.2) can be written as 



ds e - is ( E ( u )- a + XK t.t+s( u )) e - es 
To estimate it, we split K + into disjoint intervals 

h ■= [K, (i + i)0 , jeN 



(5.3) 



(5.4) 
(5.5) 



HOMOGENOUS FERMI GAS IN A RANDOM MEDIUM 



19 



of length 



C := 



\E(u)-a\ 



(5.6) 



We find 



(5.4) = 



-is(E(u)-a+\K tt+s (u)) e ~es 



je2N D ^ 

_|_ e -i(s+0(B(u)-a+A« tjt+a+c (u)) e -e(s+C)^ ; (5/7) 

where the second term in the bracket accounts for the integrals over Ij with j odd. 
Evidently, e~ i ^ E ^ u ' a ^ — e T47r = —1. Therefore, we get, for j fixed, 



e -ia(B(ti)-a+A« tlt+i (ti)) e - 



+ e -i(s+0(-B(«)-a+Xitt, t+s+c (u)) e -e(«+C) ^ 
,jj e -M(^W-«) e -e(«+() L-ttKt,t+»W _ e - iA ( s +C)Kt.t+ s (ii) ~J 

Clearly, 



^ |(5. 8)| < j dse- es e( 

j£2N + 



and 



^ |(5.9l| < / dse- £S AC 

P2Nn 



j62N 

On the other hand, we observe that for s± < S2, 



1 1 



\E{u)-a\ 



A 7T 

e |S(u) - a| 



t+s 2 



Kt )t+S2 («) - K t ,t+ Sl (u) = ( )/ ds' K s >(u) 

«t+s 2 l-t+si 

) ds' k s > (u) . 

31 ^Jt Jt ' 



1 



Since |k s '(m)| < Cq uniformly in s', we immediately obtain 

S2 - Si 



\K ttt+S2 (u) - Kt,t+ Sl {u)\ < C 



Si 



so that in particular, 



\Kt, t +c+si u ) ~ «t,t+a(w)| < C 



(5.8) 
(5.9) 



(5.11) 
(5.12) 

(5.13) 
(5.14) 

(5.15) 

(5.16) 
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Thus, we conclude that 



(5.101 < C [ ds\( 

■f- HRT ** 



(s + 



C 



This proves that for \E(u) — a\ > 0, 



(5.7(1 < 



A 



e \E(u)-a\ ' 



C 

\E(u)-a\ 



under the assumption that A = O(e). 

If \E(u) — a\ < e, then the trivial bound 



(5.7(1 < / dse'" < 



C 



is better, which ignores phase cancellations, so that in conclusion, 

C 



(5.7)| < 



\E(u)-a\+e 



as claimed. 



(5.17) 



(5.18) 



(5.19) 

(5.20) 
□ 



We may now prove Lemma |5.1| 



Proof. We consider 



~ n 

/ ds ■ ■ ■ ds n TT e 

V 1 to 



where 



and 



$(s,tt,i) := s{E{u) - a + \K t t+s (u)) 



tj := Sq + ■ ■ ■ + sj , £_x := . 
Then, subdividing K + into intervals It as before, 



i ,— ,i n e2N f ° 



(5.21) 



(5.22) 



(5.23) 



e — i*(so,«o,t-l)g-es3 _ g-i*(so+C,Wo,t-i)g— e(so+C) 



{ 



e -i$(s„,« n ,i„_i)g- 



e -i*(s„+C,n„,t„_ 1 ) e - e ( S „+C) 



(5.24) 
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we may use L — L°° bounds in Sj to get 

n „ 

( [5724] ) < H II sup / ds i 



l ,...J n £2N j=0- 



-i*(sj,«3,tj_i) — esj 



-i*(af+C.«j.*j-i) e -e(«j+0 



1 



c 



e 



" J |^(«i)-a| + e 



(5.25) 



by application of Lemma 5.2 recalling the assumption that A < O(e). In particular, 



the fact has been proven here that the bound proven in Lemma |5 . 2| holds uniformly 
with respect to t, □ 



Moreover, we prove the following stationary phase estimate. 

Lemma 5.3. Assume that f,ve iJ 3 / 2+<T (T 3 ) for some a > 0, and that A < 0(r] 2 ) 
Then, for < s < t = 0(j]~ 2 ), and uniformly in t G R. 



T 3 



due -is(E(u)- a -ie) e i\f; + ° Ks (u) ^ 



< C(a) (s)- 3 / 2 \\v\\ m/2+a{T3) || / || ff 3 /2+ „ (T3) (5.26) 



where k s = fi s * v. 



Proof. Let 



-is{E{u) — a — «e) 



T 3 



:= ' 

& a ( u ) : = e ^/ t t+s d s 'K 3 ,(«) _ ( 5 27 ) 

Clearly, the left hand side of ( |5.26 1 satisfies 

<lag s (u) h s (u) | < ||<x}- 3/2 - ct <? 5 v ||f( Z 3) || (*) 3 / 2+ff (/ft.) v ||* (z ») 
C" 

< Hfl^lU""^ 3 ) II As lliT3/2+^( T 3) (5.28) 

c 

< — ( s )~ 3/2 II / ll// 3 / 2 +"(T3)ll ^« llff3/ 2 +-(T 3 ) > (5.29) 

where we used the fact that H a (T d ), with a > f , is a Banach algebra. 
With As < 0(1), 

||/l s || H 3/2+crj- T 3-j < C ||K s || ff 3/2 + CT(T 3) < C 1 1 jU s 1 1 £1 ( T 3 ) 1 1 U | j H 3 / 2 + " (T 3 ) (5.30) 



where ||At s ||i 1 (T 3 ) < ||a*s||l°°(t 3 ) < c uniformly in s. Accordingly, (5.26) follows. □ 
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5.1. Feynman diagrams contributing to the error term. In this section, we 
combine Lemmata |5.1| and |5.2| and the results of Section |6.5| below, to prove that 
the amplitudes of all Feynman graphs that include a crossing or nesting subdiagram 
contribute only to a small error. In addition, we show that all terms involving the 
Duhamel remainder term likewise contribute to a small error. 



We point out that in our approach, the estimates provided in Lemmata 5.1 



and 5.2 and in Section 6.5 require first integrating out all time variables Sj, and 



subsequently integrating out the momenta Uj in (4.42) and (4.371 



This makes it possible to straightforwardly adopt estimates on crossing, nesting, 
and remainder terms from [5J QUI EZ] • Accordingly, our discussion in this section 
can be kept short. 

5.1.1. Crossing and nesting diagrams. We defer the discussion of nesting diagrams 
to Section [6. 5| below because the necessary ingredients are introduced only in later 
sections. We shall here anticipate the result from the analysis given there. 



Using Lemmata 5.1 and 5.2 we may infer from the analysis presented in 8, 9, 
1101117] that the Feynman amplitude of any graph n 6 L„ ^ containing a crossing or 
nesting subdiagram is bounded by 



lim \Amp w (f,g;e;<r))\ 

L— too 



(5.31) 



< || ||0 Ha || J||ooe 1/5 (log-) 4 (^ : 



V 1 



log ■ 



where In = n + n € 2N. The thermodynamic limit L — > oo and the upper bound 



in (5.31) are obtained in the same manner as in [5J O QUI 02] where we refer for 
details. We do not repeat the discussion here, and instead refer to those references. 



Let 



pc— n 
n.n 



denote the subset of Feynman graphs of crossing or nesting type, and 



Tic-n 
1 2n 



with cardinality \T c 2 ^ n \ < 2 n n\. 



U 



pc— n 



(5.32) 
(5.33) 



n+n=2fi 



Accordingly, summing over all graphs with a crossing or nesting subdiagram, 

J2 E r lim \ Am PM,9;z\ri)\ (5.34) 

* — 4 * — 4 L— >-oo 
l<Tl<N^ eT ^ 

< ll/Ha Iblla II J Woo e 1/5 ( log 1 ) 4 ( c^e'^log - ) N , 

e e 

where we replaced the factorials by AM < N . Since /, g are of Schwartz class, 

11/112, N| 2 <c. 

Moreover, 



< 1, 



(5.35) 
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from (2.22). As a concrete example, 

< J(p) = (l + e^^fP)-^)" 1 < 1, 
for the Gibbs state of a free Fermi gas, at inverse temperature < f3 < oo. 



(5.36) 



5.1.2. Remainder term. Using Lemmata 5.1 and 5.2 and following [TO] (which uses 
results in [8, 9, [17]), we straightforwardly find the following bounds on the contri- 
butions of the Duhamel remainder term. 

If at least one of the indices n, n equals N + 1, one obtains 

Jn,n), 



Km \E\p^(f,g)]\ 

L—>oo 

< ll/llallsllalUI 



r N 2 k 2 



(JV!)i/2 



max{ce-V, (ce~V) } (5-37) 



+ (N 2 K 2 e 1 / 5 + e- 2 K - N )(log-) 4 (c7 ] 2 e- 1 Nlo g 



(again using (4N)\ < (47V) 4Ar ) where the constant 1 <C K <C t remains to be 
optimized. The first term on the right hand side of (5.37) bounds the contribution 



from all basic ladder diagrams contained in the Duhamel expanded remainder term. 
A detailed discussion of an analogous result, and more details, can be found in 

on us nana. 



5.1.3. Choice of parameters and error bounds. The kinetic scaling limit asserted in 
Theorem |3.1| is determined by 



1 T 

1 = - = S 



(5.38) 



where t denotes the microscopic, and T the macroscopic time variable. Similarly 
as in [81 HOI QI|, we choose 



N = N(e) = 



log \ 



10 log log \ 
1, 



(5.39) 



K = (log^) 15 , 

so that 

e-i/ii < m < e -i/io 

k n > e- 3 / 2 . (5.40) 
One can easily verify that, accordingly, for any choice of T > finite and fixed, 

(5~34| < r, 1/15 



(5.371 < n 



1/4 



(5.41) 
(5.42) 



for r] sufficiently small. 

In conclusion, we have proven the following overall bound on the amplitudes of 



all Feynman graphs that we attribute to the error term in the expansion (4.48) 
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Proposition 5.4. For the choice of parameters (5.39), 

+ |5.37| < rj^ (5.43) 

holds, for any choice ofT > finite and fixed. The l.h.s. of (5.43) contains the sum 
of Feynman amplitudes associated to all diagrams containing crossing and/or nest- 
ing subgraphs. It also contains the sum of all terms that depend on the remainder 
term of the Duhamel expansion. 
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6. Proof of Theorem I3JJ II. Resummation of main terms 



In this section, we discuss the main terms in the expression (4.48), which are 



associated to the set of decorated ladder diagrams. This is the complement of the 
set of Feynman graphs containing crossing and/or nesting subgraphs studied in the 
previous section. We prove that the sum of amplitudes of a ll decorated ladders 
converges to a solution of the linear Boltzmann equation ( 3.4 ) , for A < 0{rj 2 ). We 



remark that because of the presence of the nonlinear self-interaction of the fermion 
held, our analysis here is much more involved than the analogous part in [51 117) . 



We let r v _ ' C r„ n denote the subset of all decorated ladders based onn + n 
vertices. Our goal is to prove that in the given kinetic scaling limit, the sum of 
amplitudes associated to decorated ladder graphs converges to a solution of the 



linear Boltzmann equation (3.4) 



6.1. Outline of the proof. Our discussion is organized as follows: 



Step 1. 



In Sections 6.2 and 6.3 we consider the Feynman amplitudes associated 
to propagator lines decorated with an arbitrary number of immediate rec- 
ollisions. The sum of all such terms produces a renormalized propagator 
(renormalization in one-loop approximation, according to standard termi- 
nology in quantum held theory). We prove that the leading term corre- 
sponds to a multiplicative renormalization of the basic propagator Ut a .t b (m)> 



see (|4.15|) for its definition. 

we discuss the Feynman amplitudes of basic ladder diagrams 



Step 2. In Section 



6.4 



(that is, the propagators between rungs do not contain any immediate recol- 
lisions). We prove that, in the kinetic scaling limit, the Feynman amplitudes 
of these graphs determine a transport equation obtained from omitting the 
loss term on the rhs of the linear Boltzmann equation (3.4). To prove that 



the error terms thus obtained are small, we apply various stationary phase 
estimates to the propagators U tay t b {u)- 
Step 3. In Section 6.6 we combine the results obtained in the previous two steps, 
and show that by decorating the basic ladder diagrams with immediate 

obtained in the kinetic 



recollisions, the full Boltzmann equation (3.4) 
scaling limit. 



As a key result of Step 2, we obtain that the dominant part of the Feynman am- 
plitudes is in fact independent of A in the regime A = 0(r] 2 ). This is a consequence 
of cancellations due to the translation invariance of the system. 



6.2. Immediate recollisions. According to step 1 outlined above, we study the 
amplitude of a progression of m immediate recollisions. 



Lemma 6.1. For any connected interval I C R+ of length \I\ < 0(t), 

du [ d se -^ t * +s ^'(B(»)-«-^ s ,(«)) 
J i 

- f du ( ' dse- ls ^ E{u) - a - K) <CX\I\ 1/2 . 

JT3 J I 



(6.1) 
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ds 



du e~ l K + " ds ' ( £ («)- A «v («)) 



T 3 



< C{1 + X\I\ 1/2 ) (6.2) 



uniformly in t > 0. 



Proof. To prove (6.1), we write 



I JT 3 



(6.3) 



(I) 



ds / due- ls{E{u} ~ a - ie} 
i Jt 3 



(6.4) 



where 

and 
(II) 

.In .It ./T3 .It 

(6.5) 

from differentiating and integrating with respect to the parameter A. To estimate 
(II) , we introduce a smooth partition of unity, 



/ / dX' ( ds ( due- is{E{u) - a - l,i) [ + ds 1 K s ,(u)e iX ' ds " K »"W , 



/ du = y2 du x 3 ( u ) 

Jt 3 ~[ Jt 3 



(6.6) 



where the smooth, T 3 -periodic bump function \j is centered around the j-th critical 
point of E(u), which is a real analytic, perfect Morse function. The properties (2.18 ) 



satisfied by v also hold for n s — v * /i s (see (5.30)). Lemma 5.3 applied on the 
support of each Xj , implies that 



t+s 



due -is(E(u)- a -ie) I dg l Ks ,( u ) e «A' , f t t+s ' ds" k.„{u) 



< C( S )- 3 / 2 e~ e 

< C(s)~^ 2 se 



t 

t+s 



ds' Xj (u) k s > (it) e l 



A' fl +s ' <Is"k s „(u) 



#3/2 + <x( T 3) 



(6.7) 



for any finite a > 0, as long as A < 0(e), and s < t 1 (where (s) = (1 + s 2 ) 1 / 2 ). It 
thus follows that 



I (II) | < A J ds (s)-3/ 2 se~ es = 0( A \I\ 1/2 ) 



which proves (6.1). In the same manner, 



T 3 



< C(s)- 3 / 2 (l + As) : 



implies (6.2). 



(6.8) 

(6.9) 
□ 
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6.3. Renormalized propagators. Next, we resum progressions of immediate rec- 
ollisions of arbitrary length, leading to a propagator renormalization. 

Let m e N, and S := {1, ... ,2m + 1}. Moreover, let < t a < t b < t. We 
consider the integral 



Tf( m ) 
U t a ,t, 



(u) := {-if) m \ ds 1 ■ ■ ■ ds 2m +i S(tb-t a -S~" s 3 ) 

« m 

/ du 2 - ■■ du 2m+ i TT S ( u 2j+ i - u^j-i ) 



'(*-,-- i-i 

2m+l 4 

TT e - '^- 1 ds '( B ^)~ A ' t s'(«^)) (6 10) 

corresponding to the Feynman graph given by a progression of m immediate recol- 
lisions. Integrating out the product of m delta distributions, 



J 2 ^dsi--- ds 2m +i Sjtb — ta — y^^j) (6-11) 

TT f du' e~ l ^- lds ' (Eiu ' ) ^ E{u) ~ HKs ' (u ' ) ^ Ks ' iu))) 
jeSn2N^ T3 

where 

tj = t a + Sl -\ \-Sj, (6.12) 

for j — 1, . . . , 2m + 1, and where 

U^ tb (u) = e -i&ds>(E(u)-x K Au)) ( 6 .i3) 

by definition of U^ b (u). 

In the following lemma, we determine the Feynman amplitude of a progression 
of m immediate recollisions, and extract the dominant part. We identify the latter 
as a multiplicative renormalization of the free evolution term (6.13|. 

Lemma 6.2. Assume that A < 0(r/ 2 ) and < t a < tf, < t. Then, for all m E N 
and v = 0(r] 2 ), 



where 

j r(main;m) 
U t a .t h 



U£i(u) = U^r m \u) + AU^l(u) (6.14) 



(«) —, (-V 2 (t b -t a ) I du> 1 . ) m ul% b {u) (6.15) 

m! V J h(w ) — hi(u) — iv J 



and 

T (m) 



\\AU£1\\ L ~ (T3) < {rnXv^ + v^){crfv^) m . (6.16) 
In particular, AZ7 t t (u) — in the case m = 0. 
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Proof. We represent the delta distribution in (6.11) by 

d 7e -<(*b-*a-E«ih 



2?r 



= e »{t b -t a ) J_ /" d7e -i(t6-fa)7 e i(E^)(7+^) ) ( 6 . 17 ) 

where we will pick f=e = j = ^rSO that (i?, — t a )v < vt = 1 and in particular, 
R v(t b -t a ) <- e un jf orm ly j n ^. j n order to keep track of the origin of these two 
parameters, we will continue to notationally distinguish v and e. 



From (6.11) and (6.171, 



U {m) (u) = e -if%ds'(E(u)-\K s , {u)+w) 

1 
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d 7e -*T(*»-*.)(_^)m f ^ . . . ds TT 

• /K + feSn2N-i 

TT f e -i!l 3 1 _ l ds'(E(u')-E(u)~ 1 -iv-\{K 3 ,(u')-K s ,{u,))) ^ 

where according to the above definitions, tj = tj-i + Sj. Our goal is to prove that 



ill 



dr) e 



-i~l(tb-t a ) 



t l \m+l / 

V 7 + iv J \ 



7 + 



du 1 



-irj 



E{u') - E(u) -j-iv 
+ erri | 



; \:\ d$'(E( u )-\ Ks ,(u)) | p y(t b -t a ) 



1 

2^ 



v 7 + iv J V 



7 + ' 



i/ t * a b ds'(B(u)-A«;,/(u)) f fa ~ ^a) r 

I to! 



du' 



-%r\ 



E(u') - - iv 
+ erri + err 2 | 

^S dU, E{u>)-E { u) 

+ cri'i + crr 2 | , (6.19) 



where we claim that 



and 



| err! | < m\v 1 / 2 (crj 2 v 1 )' 



| err 2 | < v ' (cry 
Clearly, this implies the asserted bound. 



V2 C,„2 , -l\m 



(6.20) 
(6.21) 



We first prove (6.20). To this end, we will use that || ac s || ^3/2+ (t 3 ) < c uniformly 
in s, see (5.30), and the fact that the kinetic energy E{u) is a real analytic, perfect 
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Morse function on T 3 . A stationary phase estimate similar to the one applied in 
the proof of Lemma |5.3| yields 

f du' e"*^- 1 ds '( £ ( u ') _B (")~^~ ii/ ~ A (' < v( il ')-K s '(«))) 
It 3 

du' e ~ i(tj ~* J - l)( - E( "' ) ~' E(n) ~ 7 ~ ,;,y) 



/ d\' 



< 



du' N ' ds" /t s »(0) 



—i f t J ds r (E(u') — E(u) — 7— iv— A'(k s / (u') — k s / (it))) 
6 J 



CAs, (s,-) _3/2 . 



where = + Sj. Thus, 



I er ri I < Cif m / ds t ■ ■ ■ ds 2m+1 S(t b - t a - Sj ) 



e ^>- i/2 n 



-3/2 



< Cm\is- 1/2 (Ci^V)™ . 



This implies (6.20|. 



To estimate err2, we use 
1 
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err 2 



\-y\<v J \l\>v 

du' 



(f +1 ) dl e-^'^{^-) 

( 



m+l 



- J57(u) - 7 - w 



i?7 



2m 



1 \ m+l 



\l\>v 



*r(rr-) 



2 sup sup 

d7 (HT^) m+1 



(6.22) 



(6.23) 



E(u') - E{u) -~/-iv 



- x \i\ 

' m 



(6.24) 



where 



AM 



3=1 W 



sup sup 

|-y|<2/ B(«)e[-6,6] 



1 



£(u') - £(u) - ii> 



du' 



1 

E(u') - - 7 - iv 
du' ' 



(6.25) 



E(u') - -B(u) - ii/ 



We prove below that 



A(v) < Cv 1/2 



(6.26) 
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and 



sup sup 

jGM.E(u)£[-6,6] 



du' 



E(u') - E(u) - iv 



< C. 



( 6^24) < {Ciy- 1/2 f] 2 ) m + v 1/2 {Cv- 1 r 1 2 ) m . 
Choosing v = e = 0(rj 2 ), we arrive at the asserted bound for |err2 1 



Recalling that YJj=i (™) = 2 '" - h this implies that 



(6.27) 



(6.28) 



To prove (6.26), we use that, from a stationary phase estimate, 

1 



du' 



E(u') - E{u) - 7 - iv 

f W du 'U«')-i* 



1 



E(u') - E(u) - ii/ 



u) — 7 — ii/ 

/•OO /"OO /■ 







(I 



si 



< C| 7 | / dsi / d S2 ( S2 )- 3 / 2 e-^ 



si 



(6.29) 

since I7I < v, and since E(-) is a real analytic, perfect Morse function on T 3 , as 



noted before. This implies (6.26) 
On the other hand, 

' du' 



E(u') - E(u) - 7 - iv 

ds f du' e -i<E^')-E(u)-y-iu) 



< C / ds(s] 

Jo 

< C 



-3/2 e -us 



for all 7 € E, and in particular for | 7 | > v. This implies (6.27) 



(6.30) 
□ 



6.4. Sum of basic ladders. Following step 2 described at the beginning of Section 
[6j we now determine the kinetic scaling limit of the sum of all Feynman amplitudes 
associated to basic ladder graphs. In combination with the propagator renormal- 
ization addressed in the previous section, this will allow us to complete the proof 
of the Boltzmann limit asserted in Theorem 13.11 



The Feynman amplitude of a single basic ladder graph with q rungs is given by 



wf~) (J ;/ )9 ) := { _ lf)q / dtq 



1 



*2 j-t 

dt\ I dt q 
Jo 



dt 1 



du ■■■ du q J(u q ) f(u ) g(u ) 



(6.31) 



C$K) uQ(u q ) • • • D&(«i) ^SkM <M U^(uo) 
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where the definition of U^°_ t tj (u) is given in (6.13|. 

The following intermediate result will be important for the derivation of the full 
transport equations in the next section. 

Proposition 6.3. Assume that A = 0(77), and M(ry) € N with M(rj) < 0( . '°. S " 1 ) 
and lim^o M (rj) — 00. Then, for any fixed, finite T > 0, 



£ 

q=0 



F^ c \j-f,g) := lim £ 4%T ;<?) (J; /, 5) (6-32) 



exists, and 

^(basic) ( jr \ / j ~T7 \ / \ rp{basic 



Fr ' (J;f,9) = J du f(u) g(u) Fr SlC ' (u) (6.33) 
where F^ aslc \u) satisfies 

dT p(fiasrc) {u) = j ^ ^ _ ^ ) ^asic) £ u ,j _ (g 34) 

wii/i initial condition Fq(u) — J(u). 

This proposition is an immediate consequence of the following lemma. 
Lemma 6.4. Assume that X = o(rj). Then, 

7 ,{basic;q) / T r \ 7 Abasic— main;q) , T r \ 1 r-\ 

u\ '{J;f,9) - u\ {J;f,g) (6.35) 

< CqXt 1 ' 2 {rftlog-y- 1 + Cqfji^tlog-)"- 1 , 
rj rj 

where l4f baslc mam ' q \j] f,g) is defined in \ 6.J^ r l) below, and 

F (ba^ q){J , Lg) . = UmW (^c-mam; g ) (J;/)ff) ( g 3^ 

exists for any fixed, finite T > 0. In particular, 

F<r basic -' q) (J;f,g) = / du/MsW^^W (6-37) 



T 

where p^ aslc ><f> ( u ) satisfies 



and F^ baslc ' q) (u) =0ifq>l,and F (baslc ' 0) (u) = J(u). 



Proof. First, we write 



so that 

^ 1 ,*,fe)^ 2 , t ,_ 1 (%--i) 
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and 



r(0) 



(0) 



r(0) 



where £/ t ^(u) = ul?\(u) (due to J* t * = — J.' on the rhs of (6.13)) independently of 
t > t' or t < t! . Accordingly, we find 



LI 



(basic;q) 



(J;f,g) 



(-rfy / dt q / dt q ^ 



dtl / dt q / 



«2 



dh 



du ■■ ■ du q J (uq) f(u ) g(u Q ) 



(6.42) 



where, evidently, {u q ) = 1. Let 5„(u') := 5(u' — u) so that 



(basic;q) 



(J;f,g) = / d W J( W )W, 



,(basic;q) 



(6.43) 



Our proof comprises the following main steps. 



Step (1) First, we verify that u[ haslc,q \j; f ', g) satisfies the following approximate 
recursive identity, 

= -rf du q J(Uq) dUq^i 



dt„ I dt 



+ o( t 



,-i/ 2 {Cr, 2 ty 

(9-2)! 



(uf t {Uq)U?\(u q ^) 

(6.44) 



We note that the main term in (6.44) differs from (6.42) only by the upper 
integration boundary for the variable t q —i, which is replaced by t q _\ — > t q . 

Step (2) Next, we prove that for A = 0(77), the nonlinear self-interaction of the 
fermion field only contributes to a small error, 



LI 



(basic;q) 



{J; La) 



rj 2 du q J(u q ) diiq-i / dt q I ds exp ^ - is (E(u q ) - E(u q -i))^ 



+ 0{Xt 1 ^( V 2 t\og-) q - 1 ) + 0(r ] ( V 2 tlog-) q - 1 ). 

77 V 



(6.45) 
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Step (3) Iterating (6.45), one gets 



7 .(basic-q) / T r \ 7 / ( basic— main:q) / T r \ 

u\ {J; f,g) = K J {J]f,g) 

+ 0(qX < 1/2 ( rft log - Y- 1 ) + 0( q r] ( rft log - ) (6.46) 
77 77 



where 



7 .(basic— main-q) / T r \ 

Ul J {J;f,9) 

'■= (-V 2 ) q J du q J(u q ) J du q _ 1 ■■■ du a f(u ) g(u ) (6.47) 

rt />t 2 rt pt 2 Q_ , 

j dt q ■■■ I dt\ I ds q ■■■ / dsi TT cxp ( — isj(E(uj) — E(iij-i)) 
Jo Jo J-t J-t 2 =1 ^ 



Introducing the variables T = r] 2 t and Tj = r/ 2 tj, 

7 .(basic— main:q) / T r \ 

= (-1) 9 j du q J(u q ) diiq-x ■ ■ ■ du f(u ) g(u ) / dT q ■ ■ ■ j dT x 



j-T/n 2 rt 2 /r, 2 1 

\ ds q --- / dsi TTexp ( - isj(E(uj) - E(uj-i)) ) . (6.48) 

J-T/n 2 Jt 2 /;f- X = l V ' 



For every T J+ i > 0, 

T 3 + i/v 2 . . 

d Sj exp -is^^J-BM) — > <5(£?(u J )-S(u i _ 1 )) (6.49) 

weakly in the limit 77 — > 0. Therefore, 
lim^™ in;9) (J;/ 5 5) (6.50) 



(-T) 
9! 



r 9 

/ dn ••■ du q J(u q ) f(u )g(u ) J| 5(E(uj) - E(uj_i) ) 
^ i=i 



It is easy to check that (6.501 has the properties asserted in the lemma. For 



more details on the limit 77 —> 0, we refer to [IT]. 



It thus remains to prove steps (1) and (2). 
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Proof of step (1). The difference between uj: baslc ' q \j; /, g) and the first term on 
the rhs of (6.44) is given by 



MA, 



(basic;q) 



( First term on rhs of §6A4§ ) - u[ haslcvl) (J; /, g) 



(~V 2 ) q / dt 



dt 



9-1 ' ■ 



dt x I dt. 



q I dtq-l ■ ■ ■ 



1-2 



dtt 



du ■■■ du q J {u q ) f(u ) g(u ) 



(6.51) 



The only difference between this expression and the expression (6.42) consists of 
the integration boundaries for the variable t q -i. 



We use a stationary phase argument similarly as in Lemma |5.3| to bound the 
integral in u q , which yields 



dUqJ(Uq)Uj°l(Uq) < C || J || fff + „ (tq - tq) 



(6.52) 



for some a > 0, and similarly, 



du f(u )g(u )U^ ti (u )\ < Cll/H^ || fl || K§+ . 



recalling that Hi +a (T 3 ) is an algebra. 

For the integrals in Uj with j = 2, . . . , q — 1, 

r(°) c.aTTW / 7. c+. . 7. .w-3/2 



Accordingly, writing 

B J,f,g II ^ ll H |+- (T 3) II / ll ff |+- (T 3 ) II 9 ll H |+- (T3) 



(6.54) 



(6.55) 



we find 



AU (^c; q){J . Lg) 

< C^Bjj^r, 2 " f dt q f 
Jo Jo 



dt„ 



1 



d<„_i / dt 



(t q ~t q ) Z/2 
1 



(t, - *, - (t q -l ~ t 9 -l)) 3 / 2 



dtq-2 



to 



dh 



9-1 

n 



tq-l 



dtq-2 



dt 1 



(tl-tl^fJiiti-tj-fa-l-ti-l)) 3 '*' 



(6.56) 
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To bound the integrals in tj with j = 1, . . . , q — 1, we use 



tq-l 



(It 



9-2 



(t q -l - t,_l - (t,_ 2 - V2)) 3/2 



*2 



r/fi 



1 



< 



rff 



(t 2 _ *a - ( tl - fl ))3/2 (ti-t^a/a 
1 



9-2 



(Vi-Vi - (i,_ 2 -t,_ 2 )) 3 / a 

1 1 



f/fl 



< 



(t 2 -t2-(ti-ii))8/ 2 (tl-il) 3 / 2 

1 



( / ~„\/n V 3 SU P / ^1 ' 



(*i -*') 3/2 (*i} 3/2 



< C q ~ 2 , 



(6.57) 



where we first majorized the expression by extending all integration intervals to 
K, and subsequently translated all variables tj — > tj + tj. The integrals in tj with 
j = 1, . . . , q — 2 are easily seen to be bounded by / q _ 2 \\ ■ We thus find 



Awf~>(J;/, 5 ) 
< C« Sj./.a ?7 4 / di 9 sup f f dt q 1 

JO tg JO \tq — t 



(t q -t q ) 3 / 2 



dtg-i SUp 



dt q -l 



4 



"^{t q -t q -(t q . 1 -t q . 1 ))^ 



(9-2)! 



i 2 *19-2 



; (r ? 2 *)9-2 



(9-2)! ' 



using that 



dt q -i sup 



dt„-i 



(t q -t q - (t,_! - Vi)) 3/2 



< 



Consequently, 



AW (^c; 9)(J . /)5) 



-1/2 



(9-2)! 



(6.58) 



(t q -t q ). (6.59) 



(6.60) 



follows, as claimed. 



Proof of step (2). Similarly as in Lemma 6.1 



( U^ tq (u q ) («,_!) ) = exp ( - i{t q - t q ) ( E{u q ) - E(u q ^) ) ) 

x ( 1 + g(u q ,u q -i;t q ,t q ; A) ) , 
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for g{u q , u q -i,t qi t q ; A) = e *« - 1 satisfying 

ll<?( • ; u q~l\tq,tq'i ^)\\h 3 / 2 +° < Co A I ^ ~ *g I 



(6.62) 



as a function of u g , with Co dependent on a > 0, but independent of i g , i g , A. 
The regularity of g is inherited from the pair interaction potential, v <G H 3 ^ 2+a (T 3 ) 
(see (2.18)), and proven as in (5.30). 



Similarly as in the proof of Lemma |6.1[ a stationary phase argument yields 
du q exp ( - i(t q — t q ) (E(u q ) - E(u q -i) ) ) g(u q ,u q -v, t q -i, t q ; A) 

< CXitg-tg)- 1 / 2 . (6.63) 
For the integral in t qi we thus find 



dt„ 



du, 



q (uf \ (u q )uf\ (« g _o 

V L q i L q L q j^q 

< CXt 1 ' 2 . 



It is then straightforward to arrive at (6.45). 



(6.64) 
□ 



We may now prove Proposition |6.3 



Proof. To establish Proposition 6.3 it suffices to verify that 

M(ri) M(n) 



E U$Jp q \j;f,g) Y, = 0. (6.65) 

g=0 g=0 

For the left hand side, we obtain the bound 

M(n) 

£ |4/? c;9 V;/,ff) -4/" Qm;9) (J;/,ff) 

M(tj) 

< 

"77' ' 77 



£ [At 1 / 2 ( CrjH log i) 9 " 1 + ^(C^logi) 9 " 1 

9 =1 



(6.66) 



using the estimates on the error terms in (6.46). For the choice of parameters 

(6.67) 



A = 0( V 2 ) , t = 0(?T 2 ) , M(rj) 



where we assume that cq > 2, this is bounded by 



logi 



Cq log log i 



(6~66l) < 7] (C log - ) M ^ M(ri) 
i r v 



< i] ■ r) 

i 

< ?y 10 



, log Op . i \ -. 

" ( loglog I , 1 

1 log - 
»7 



if 77 is sufficiently small. This implies (6.65). 



(6.68) 
□ 
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6.5. Remark on bounds related to nested diagrams. The analysis given 
above enables us to control the Feynman amplitudes associated to nested diagrams, 
as was noted in Section |5.1| To this end, we recall that for the Feynman ampli- 
tudes belonging to all diagrams containing crossings or nestings are evaluated by 
first integrating over all time variables Sj, and subsequently integrating over mo- 
mentum variables Uj . We note that our argument below exhibits the same ordering 
of integration steps. 

Following the definition of nested diagrams in Section |4.2| a nesting subgraph 
of length m is a progression of m consecutive immediate recollisions connected via 
a propagator line to two outermost vertices which are mutually contracted. The 
contribution to the Feynman amplitude associated to this segment of the graph is 
proportional to 



T 3 



du I ds U, 



(6.69) 



for some r < 1 = T/rj 2 , where the integration variable u appears nowhere else in 
the full expression of the Feynman amplitude (for the entire graph) . The factor rj 2 
accounts for the contraction of the two outermost vertices. 



Using Lemma |6.2| one straightforwardly verifies that 



rf I du 

/T 3 



< 



dsU 



(m) 



t a ,t a + s( U ) 

v 2 [ du fdsU^?\u) +v 2 tsu P \Au£1 +s (u)\ 
Jt 3 Jo s , u 



< C m TJ 



ro ~ 2 ' du 

T 3 



+ c m V: 



(6.70) 



where we have recalled that in Lemma |6.2[ the values of the parameters are given 
by A,^ 0(77 2 ). 

We claim that 



\E(u)\ 



"1 rnj 



(6.71) 



The proof of (6.71 ) is similar to the one of Lemma 5.1 If |-E(w)| < r 1 , the bound 



is trivially fulfilled. If |-E(it)| > r 1 , we define 
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and the intervals Ij := [j'C , (j + 1)C) with j € J := N n [0, J]. Then, the left hand 
side of (6.71 1 can be bounded by 



U, 



(0) 

t OJ t a + S 



(«) 



(a + C)' 



(0) 

t Q ,t a + s 



(«) 



je2N n./ , ' J i 



i ta+jC,t'+jC + »( 1i ) 



V / ds e - isE{u) ( e isXKt 

- fs 4- f*V 

_ e *(»+0^"t o +t f +i)c,t'+w+i)c+*C u ) v ' ^ 



j'e2N nJ 



Similarly as in the proof of Lemma |5.1 



3 is\K ta+jCtta+j(+3 (u) S e l O+C)^t a + (j + lK,t Q + o + lK+ s (")| < (7AC- 



and evidently, 



Accordingly, 



< 



(s+cr 



< c 



C 



(a + C) m ~ 
(m- 1)! 



7o (m-l)! 



(2r) 



for A,r _1 < 0(r) 2 ). This proves ( |6.71 ). 
Therefore, we conclude that 

(0.70) ; C"V / </</ 
1 



E(u)\ + V 2 



< C m r, 2 log - + C m r, 
V 



(6.73) 



(6.74) 
(6.75) 



(6.76) 



(6.77) 



The gain of a factor r\ is crucial, and immediately implies the bounds on nesting 
subgraphs used in Section |4.2| For a more detailed discussion of nested diagrams 
in the context of the weakly disordered Anderson model, we refer to [T71 [S] ■ 



HOMOGENOUS FERMI GAS IN A RANDOM MEDIUM 



39 



6.6. Decorated ladders and Boltzmann limit. Following the list of steps ex- 
plained at the beginning of Section [6j we now carry out step 3. Combining Lemma 
|6.2| and Proposition |6.3| we derive the Boltzmann limit for the sum of decorated 



(renormalized) ladders, and complete the proof of Theorem 3.1 



For notational convenience, we introduce the multiindiccs 



,<«) 



(m , • • • , m (9) := (m , • • • ,m q ) 



(6.78) 



for fixed q £ N, and 



|m^| := m + ••• + m q , |w | = m a 



(6.79) 



We use N(e) — 0{ l( ^°f g i ) with ~ = i = ^ as in (5.39), and consider 



: = E E E r lim Am P 7 r( J ;/)fl';*; 7 ?) 



n=0 n+n=2n wC r (ia - d) 



E E E 

M ,...,M JV(e) eN g gN m,C<!),m(a) 

E«j<"(«! | m(«) I +|m (8) |+<,=M„ 



{ren\q\m} q ' ,m ) / 7. 



(6.80) 



where 



14™**** ^ \j\f,9) ■= (~V 2 ) q I du ---du q J(u q )f(u )g(u ) (6.81) 

8 



di 9 • • • / dti / cft 9 • 

JO JQ JQ 



3=0 



using the convention that i_i = = This is the Feynman amplitude of a 

ladder with q rungs, where the two particle edges labeled by the momentum uj are 
decorated with rrij, respectively fhj, immediate recollisions. 
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To extract the dominant terms in this expression, we define 

(ren—7nain;q;rn^ q ^ ,mS q ^ ) , 



\J;f,g) 



12 



{-rj 2 ) q dt q --- dtt dt q --' dti / du ■ ■ ■ du q J(u q ) f(u ) g(u ) 



(-7] 2 ) q dt q -- dti dtq--- dti I duQ ■ ■ ■ du q J (uq) f (u ) g{u ) 



3=1 



j=i j 



TO,-! 



E(v,j) — iv 



(6.82) 



where 



U^_ ut .{u r ,u ) := U^ ut .( Uj )U^ ut .{u ), 
and where we have inserted a factor 



r(0) 



r(0) 



(6.83) 



(6.84) 



To get (6.82) from (6.81), we have replaced the contributions from the immediate 
rccollisions by their dominant parts identified in Lemma |6.2| 

Moreover, we define 

i(ren—7nain—Q;q;mS q ^ ,m^) , 



u 



'(■/;/, 5) 



{-rj 2 ) q I dtq-- I dt\ I dtq--- I dti / du ■ ■ ■ du q J(u q ) f(u Q ) g(u ) 



(v!) — E(uq) — iv 



(6.85) 



77 {tj-tj-i) / 



£(u') - -B(wo) + if 



yn 3 1 



To get this expression from (6.82), the momenta Uj have been replaced by uq in 
the last product in (6.82). Here, we anticipate the conservation of kinetic energy 
in the collision processes, which will emerge in the kinetic scaling limit. We then 
prove the following result which immediately implies Theorem |3.1| 
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Proposition 6.5. Let N(e) be as in Section 5.1.3 Then, for any fixed, finite T > 
and t — , and v = e = \ , 



E 



lim 

Mo,-,«iv(<)e B o gSN 
EMj<»(0 



E 



\ m te) \+\m (q) \+q=M q 



7 , (ren— main— 0:q:m^ q \m^ q ^) / T r \ 
T/rj 2 ~ ( J 'Ji9) 



duf{u)g(u)F T (u) (6.86) 



where Ft(u) satisfies the linear Boltzmann equation (3.4-) with initial condition 
F (u) — J(u). Moreover, 



rv^M E E E 

M ,...,M JV(e) eN ggNo m(3),„(8) 

E Afj <W(e) |m(«)| + |m(9)|+ 3 =M<, 



4/7 9! - <,), - 5)) (^/>3) 



(6.87) 



7 Aren— main:q:m^ q ^ .m' 9 ' ) / T /• \ 

-^4/^ ~ (J;f,g) 



and 



lim 



E E E 

M ,...,M JV(e) eN qgNo m(8),m(a) 

fifW | a (9)| + | S M| + J=M , 



(ren—main;q;m^ ,mS q ^ ) 
T/r, 2 



(J;/, 5 ) (6.88) 



7 /(ren— main— 0:q;m^ q ' ,m' 9 ' ) / t p \ 
T/n 2 (J;f,9) 



Proof. We first verify the Boltzmann limit for the main term before proving the 
error estimates. 



1. Proof of ^J^ . We have 



E E E 



jTo ^ ^ ^ 

M ,...,M^ (e) €N g eN jnta),mt<!) 

EM 3 <«(«) |H («) |+ |g(j)| + , =Ml 



7 Aren— main— 0:q:m^ qj .m^ q! ) / T p \ 



E E E 

Mo,...,Af;j,—eNo <?<EN 



7 /(ren— main— 0:q:m^ q ^ .mS q ^ ) / T r \ 

U y Th 2 - '(J;f,g) 



E E 



mil) 

m(«)| + |m(«)| + ,=M, 

UX,,^ ""- - '{J;f,g) 



l T/r, 2 



(6.89) 
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We find 

j(ren— main— Q;g;m ,m ) f j t 



E 7 Aren— main— 0:q:m^ q \m'' q " > ) / T /» \ 
(J',f,g) 



rt pt2 ft ft<2 

{-rj 2 ) q dti--- dtx dt q --- dt l j du ■ ■ ■ du q J{u q ) f(u ) g(u ) 
Jo Jo Jo Jo 



[] Uif_ utj ( uj ;u ) Ug\~i ur ,u ) (6.90) 

7=1 

1 _ _ 



First of all, we note that the expression obtained from setting the product on the 
last line equal to 1 is precisely U^ baslc ' q \j \ f, g)\ that is, the amplitude of a basic 
ladder with q rungs from Section |6.4| Clearly, the product on the last line equals 



e ^ 12te J fe ' gTO , ( 6 . 91 ) 
We are choosing v — e, and consider the limit v = e — > 0, where 

lim Im f du' 1 = 7T fdu' S( E(u') - E(u) ) . (6.92) 

v-*o J hj\u' ) — kj(u) — iv J 



Thus, it follows straightforwardly from this, and from Proposition 6.3 that for 
A = 0(r) 2 ), and any fixed, finite T > 0, 



j^iren—main—O)/ T r \ v \ 7 Aren— main— 0:q) / T /. \ /~ 

Ft (J;f,9) ■= 1™ 2^ U T/ni (J; 1,9) (6-93) 

17 ?eN 



exists, and 



F T r —°\j;f,g) = j duf(u)g(u)F T (u) (6.94) 



where 

F T {u) = e- 27rT / du ' KE(u)-E(u')) phasic) (u) ( g Q5) 



(see Section 6.4 for the definition of F T baslc \u)) satisfies the linear Boltzmann 



equation 

d T F T {u) = 2tt / du S(E(u) - E(u ) )(F T (u) - F T {u) ) , (6.96) 



with initial condition 

1 



F (u) = lim J(u) = lim ~ p ( a+a u ) . (6.97) 

L—tOC L—*OC ±j 



This proves (6.86) 
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2. Proof of (6.87 ). Recalling Lemma 



6.2 



we consider 



U 



T/v 2 



(J;f,g)-U^ 



(ren—rnain;q;mS q ' ,rn^ q ^ ) 



du - ■■ du q J (u q ) f(u )g{u ) / dtq- ■ ■ dti \ dt q - ■■ 

Jo Jo Jo 

3=1 



(6.98) 



dti 



(A) + (B) (6.99) 



where 



(A) := 



1 r 

(~V 2 ) q ^2 I du o---du q J (u q ) f(u ) g(u Q ) 
1=0 J 

t />t2 rt i>t2 ^ 

dtq ■ ■ ■ dti / dtq ■ ■ ■ dti 



3 = 1 



3=1+1 



(6.100) 



The functions £7, 



6.2 



t„ tT™' m ^( u ) and ^Ut™l(u) were defined in connection with Lemma 
(i?) is the analogous term with the roles of the variables and tj exchanged, 



and with 



ith ^jU-W re P laced b y u t™-™P'( u i) 



(main;m) , 



To bound the integrals with respect to tj in (6.100), we recall from (6.10) that 



(-?7 2 ) m / ds\ ■ ■ ■ ds 2m +i Sjtb-tg-^Sj) 

J R 2rn + l 

~ m 

/ du 2 ---du 2m+ i T\S(u 2 j+i -u 2 j-\) 

J (T 3 ) 2m J"_, 



2m+l 

n 



"*'/*/_ 1 rfs'(-E(u;)-AK s , 



(6.101) 
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Thus, 



1 rU-U-i 



jdu a ---du q J Q dt q ---J o 2 dt, J] ( (uj) U£l~ (u Q ) ) 

duo ■ ■ ■ du q J j 

n " -j 

IT (-V 2 )™ 3 / dUj j2 • • • f/» ,.2/-. • I TT <5(uj,2i+l - Uj,2i-1 ) 

7=1 • / ( T3 ) 2m3 f=i 

/ 2m .+i ' ' ' <%2m 3 +l - XI S 3V» ) 



+ 

2m 3 +l 

n 



,-*/«/_! d*'((B(«j,/)--B(«o)-A(K i /(tt j ,/)-»s i /(t»o)) 



(6.102) 



We now bound this expression by performing the integrals over all time integra- 
tions first, using Lemma |5.1[ where we begin with the latest time, and successively 
integrate out the preceding time variable. We obtain that this is bounded by 



|(6.102)| 



r q r 3 

< / du ■ ■ ■ du q TT / du j}2 ■■ -duj 2 m , + i TT S(uj t2i+1 - %, 2 i-i ) 



2m j +1 



3=1 



2\m.j 



n 



i 



\E{ Uj )-E{u )\+u 1J- |£?(^)-£(V)|+zv 



(6.103) 



Integrating out all delta distributions, we obtain 



(6.1031 



< 



du 



E(u Q )\ + v 
1 



fflj + 1 



\E(u) - E(u )\ + v 



< f[(c v 2 ^\\o g -) 2 r 3 , 



(6.104) 



3=1 



for v — e = 0(t7 2 



From |Z/A U ^ _ ( u )| = 1, and ( |6l5| ), we find 



rrijl 



so that 



II KOI < e C " < e c " 4 < C. 

3=^+1 

since T = rj 2 t is fixed, and hence of order 0(1) in rj. 



(6.105) 



(6.106) 
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Moreover, Lemma [672] implies that 



< eCv^tj-t,-!) + rjC m :l ! 
(where we may, for instance, assume c = 2c' > 1) so that 



(6.107) 



3=1 



3=i 



< e^'^Cti-tj-O JJ(i + ^c™,) 



3=1 



3=1 

< e c " 2t (l + 2ryc M ') < C . (6.108) 

For the second inequality, we used e^ 2 ^ - *^ 1 - 1 > 1 and c > 1. For the last 
inequality, we have recalled that M g < N{e) = ^ [tog^i » fr° m the beginning of 
Section [6T6| This implies that 

r\c Mq < 77 77 _Cl °sii°E" < (6.109) 

for ?7 sufficiently small. Moreover, T = rft is fixed and of order O(l) in tj. This 
implies that (6.1081 is bounded uniformly in 77, for 77 sufficiently small. 

Hence, we conclude that 

\(A)\ < (c v 2 yJ2 fto*-~ C dtiWAooWfWooWgU 

e=0 Jo Jo 

lAC/^MlfK^-^logi) 2 )'' 



(6.110) 



Using the bounds (6.16 1, 



3=1 



I (A) 1 < (a^iiJiUH/iuiifliu V fdt q --- f 2 dh 

e=o Jo Jo 

n (c^oogi) 3 )^, (6.111) 

3=1 ^ 

< (Criogij'HJlUII/IUIIfflloo ^ 

f[ (C^-^log-) 2 )^. (6.112) 
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We notice the crucial gain of the factors \v~ x l 2 and v 1 ^ 2 . 
The term (B) can be estimated in a similar way. 
In conclusion, we arrive at 

EE E [(a) + (b) 

M ,...,M JV(<) eN ggNo m(9),m(9);|m(«)|+|m(«'|+g=M„ 

EAf 3 <JV( £ )/2 — — — — 

< E E M i i Xv ~ 112 + v 1/2 ) (Crj 2 ^ 1 {log ^f) M " 

M ,...,M N(e) eN q<N(e) 

EAfj<N( E )/2 



< {2N{e) + 2)\{\ V - 1 ' 2 + v 1 ' 2 ) {Crf v-^log-) 2 )"^ 



< V 



1/10 



(6.113) 



for the choice of parameters of Section 5.1.3 and for v = e. This proves ( 6.87| ) 



• 3. Proof of ( 6. 88 ). To begin with, we note that 

1 



du' 



du'- 



E(u') - E(uj) -iv J E(u') - E{u ) - iv 

= (E(uj) - E(u ))m,jGi(uQ,Uj;v] 



(6.114) 



where 



du' 



From a stationary phase argument, 



E(u') - E(u ) ~ iv E(u') - E(uj) - iv 

d Sl ds 2 [ du' e -isi(E{-u')-E(uo)-iv) e -i S2 (E(u')-E(uj)-iu) _ 



(6.115) 



I Gi(u ,Uj-,u) I < / d Sl d S2 ( Sl + S2 )- 3/2 e-^ +S2) 

JR + xR + 

ds' ds{s)- 3 / 2 e- vs 



s>s'>0 

< [ da{s)- l '\ 

Js>0 

: Cv- 1 ' 2 . 



(6.116) 



More generally, for rrij £ N, one can straightforwardly show along the same lines 
that 



n( 



du'- 



1 



"-n(/""'; 



1 



E{u') — E(uj) — if J 7-o E(u') - E(u ) -iv 



J2(E(u t ) - E(u ))G e . mM (u^;u) 



(6.117) 



1=0 
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with := (ito, . . . , u q ), for functions 



G, ;m( „(M (9 V) := (n^wrO ( n (^(«o)) 



(6.118) 



i'=0 



with Ri U (u) := f du' E r u /\_E( u )-iu - < - >ne easu Y sees that, for I G {0, ... , 

|G £; „(„(m (9) ;^)| < m/C^-^i/- 1 / 2 , (6.119) 



for a constant C independent of £, v and {to,}, using (6.27) and (6.1161 
Next, we observe that 



(E( Ui )-E(u Q ))\j\f_ ut .(u r ,u Q ) 



= {id Sj e 

where Sj = tj — tj—\ > 0. 
Therefore, one finds 



isj(E{uj)-E{u )) j e iA /t/_! «V («;>)-«<*' («o) ) 



(6.120) 



7 Aren— main:q:m^ Q ^ .m*- 9 ' ) / t p \ 7 Aren— main— 0:q;m^ q \m^ q ^ ) / T r \ 

U\ '- '(J;f,g)-Ul - '(J;f,g) 

= (I) + (II) + {III) + (IV) (6.121) 



with 



9 ,.* 



(/) := i(-r, 2 ) q Y. / d *- 



£=0 











/o 





\ / \ rrtj 



(6.122) 



du 



1 



E(u')- E(uj)+iv 



te+i—ti- 
s e =0 
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and 



(II) := ^h 2 )'^ / dt q 



dsj 



dtj^2 / dti 







f dt q -- 




h 





r] 2 (t e - U-i) I dv! - 



m£< 



'/ J 



3-1 



i— 1 

(rj 2 Sj) m :i ^\J^_ x ds' (K B ,(uj)-n e , (u )) (6 123) 

{ ft ^l*,to;^)^(-^&-*w)r}" ■ 

where for each integration by parts has been applied to the variable S£. With 
respect to the latter, the expression (6.122) comprises the boundary terms. 

The terms (III) and (IV) are similar to (I) and (II), but in (177) and (IV), 
integration by parts is applied to the variables ~Sj = tj — tj_\. Accordingly, the 
roles of Sj and ti are exchanged with those of Sj and U, respectively, and moreover, 
W du ' E (u>)-E(u 3 )^ y nj is exchanged with UU du' EW _i (tlo) _ to ) TOj . 



Bounds on term (I). Clearly, 



|(/)| < V 2q II J||oo ll/lloo llfflloo SU P \G e , m ^(u^;^)\A 1 A 2 (6.124) 



where 











/•*2 


Ai := 








■ / dti 








Jo 


Jo 



ni 



tj — i,tj 



(6.125) 



1 



E(u') - E(uj) + iv 



and 



A2 := sup / dti ■ ■ ■ 

u ,--- ,u q Jo 



tt 



dt e+ i I dti^i ■ ■ ■ I dh 







1-2 



(I 



n{i(^-'->n 



sc=0 



(6.126) 
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We use (6.71) and bounds similar as in the case of (6.102), to estimate the factor 
(6.125) involving the integrals in tj and Uj. Thereby, we obtain an upper bound 

Ax < (dog-)S^ . (6.127) 
v 

Furthermore, we note that for every fixed index i in the sum, there are only q — 1 
integrals with respect to the variables ti , keeping in mind that there is no integration 
over tg. Accordingly, we bound the integrals in tj in A2 by ^~]Tj (the gain of a 
factor j as compared to t q is crucial). Moreover, it is evident that 

q -. 

E II ^(v 2 (tj-tj-i)) ' = = e"'* < C (6.128) 

mo,...,ro s j=l 



holds for the integrand in A2. Recalling the bound on \G t . m ( q ) {yS q ^\ u)\ in (6.119), 
we straightforwardly obtain 

E 101 

l<|m('!)| + |m< < ')|<M, 



E 

l<|m(9)| + |m(<J)|<Af, 
M 2q , 



[q-l)\ v 



/2 (Crr^\ 1 
(2«-l)!" (9-1)! ( ° l ° g J 



(6.129) 



where the factor q accounts for the sum with respect to I = 0, . . . , q, for each fixed 
m( gS> and mr 9 \ Moreover, we have used the estimate 



#{(' 



7(9) rnM 



) j \m {q) \ + |m^| < M q } < C J2 



r 2q-l 



CM 2q 



< 



r<M„ 



(2q-l)\ (2q-l)\ 



(6.130) 



from 4{{mS q \m (q) ) \ |mW | + |m (, » | = r} < C-^-^ 
of lattice points in a simplex in N^ q of side length r 



the latter bounds the number 



Hence, we conclude that, for any fixed T — r/ 2 t > and r\ sufficiently small, 

l<|m(3)| + |m (<l) |<A/ 



< qr/e 



-3/10 



rjiq- 



(9-1)! 



1/2 



(6.131) 



using that M a < N(e) = '° g , g — r (see Section 

° y V ' 10 lo£f log - v 



5.1.31, v = e = 0{if), and g > 1. 
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Bounds on term (II). For the term d Sj [ • • • 1 in (6.123), We note that d Sj f(tg) 
f'(ti) for every I > j, since te = s + • • • + Sj + ■ ■ ■ + S£. For i > j, 

_ Q g-'/tV-i da'(-B(«i)--B(«o)-A(it,/(ti«)-K i /(t*o))) 

Sj 

AO) 



= iX (« a '(ixi) - K S '(tto)) f7 ti _ i4 (ui;uo) 

S 7 — ti — 1 



(6.132) 



There is no term proportional to (E(ui) — E(uo)) on the last line because d Sj ti 
1 = d Sj ti-\. Similarly, 



(6.133) 



whenever i > j. Using the a priori bound || n s ||i/x<nr3) < c i uniformly in s, (a 
consequence of the fermion statistics, as we recall), we conclude that 



d, } [ •••] in (6.1231 < C\[ 



(6.134) 



K - 1)! 



+ A- 



(mj)! 



where < fj - ti_ x < t = 0(r/" 2 ). Combined with ( 6.102 ), we arrive at 

M 2q 1 
£ \{U)\ < ^-^(n 2 + \)(ci 1 2 tf (c\og-) M ^ 

l<|m(9)| + |m(«'|<M g 



(6.135) 



where the gain of a factor (rj 2 + A) is crucial. Using the same arguments as above 
for (16. 1311), we find 



£ K JJ )i < ^ 1/2 

l<|m(")| + |m (<j) |<Af, 



(6.136) 



for every fixed T = r/ 2 t > 0, for r/ sufficiently small, given that M q < N(e) 
ioio g S iLi ^ = e = (^)' and?>l. 



Concluding the proof. The terms (///) and (/V) are estimated similarly, and yield 
similar bounds as those derived for (I) and (//), respectively. In conclusion, we 
find that for any T = rj 2 t fixed and rj sufficiently small, 



£ E 



E 



(/) + (//) + (///) + (IV) 



JV/ ,...,Af JV(e) eN qgN m,(s) 

2«j<«W |m(9)| + |m(8)! + ,=M 9 



< (7V(e)) JV ( £ )4r / 1 /2 £ x(«<iV(e)) 

< 2 7 7 1 / 2 (7V( e )) JV W+ 2 



< V 



1/3 



(6.137) 
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that is, with A, v — 0(t] 2 ), and 



for the choice of parameters given in Section 5.1.3 

= TmT^rJ^x- Here, we have used the trivial but sufficient bound 



10 log log | ' 

#{ (Mo,-- 



) e < (c) | m j < N ( e ) } < (N(e)) N(e) , 



(6.138) 



and observed that in the sum with respect to q, the condition s M j < N(e) and 
+ q = M q imply that q < N(e). In conclusion, (6.88) follows. □ 



This completes our proof of Theorem |3.1| 
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7. Proof of Theorem [372 



Based on our p roof of Theorem |3.1[ the proof of Theorem |3. 2 1 is straightforward. 



Applying Lemma 5.1 with £ = A 1+d = r/ 2 , and <5 > arbitrary but fixed, we 



immediately conclude that the estimates for crossing and nesting diagrams used in 



the proof of Theorem 3.1 remain valid 



Adapting the proof of Theorem 3.1 to the relative scaling of parameters as as- 
serted in Theorem |3.2[ yields the bound 



(7.1) 



\\m jAmp, (/,<?;£;??) | < C(JJ,g) (logi) 4 (ci 7 2 A- 1 nlogi) fl 



L— f oo 



(see (5.31) for comparison). Since A = 77 , we find 



V V lim I Am Pl (/,j;e;))) | 

* — J — ' L — >oo 



l<n<N nGT 2 f\ 



(7.2) 



< N\(logj) 4 (cX s Nlogj) N 



In contrast to (5.34), we are not carrying out any classification of Fcynman graphs. 



The sum (7.2) extends over all classes of graphs, including ladder, crossing and 



nesting graphs. 

For the given sc aling of parameters, we note that instead of the overall factor 



e i/5 ~ ^(i+iVio in ( [5.34^ W e are now obtaining a factor r/ 5N which is < A^+^Z 10 
for A sufficiently small, and for parameters chosen similarly as in Section |5.1.3| 

In particular, we can now substitute A~5~ for 77, in the bounds given in Section 
so that N = O L to f*i ). 

v log log 4 I 



5.1.3 



Likewise, the estimates on the Duhamel remainder term of Section [5.1.2| can be 
easily adapted to the present case, and we find that 



dpf( P )g(p)F T (p) = lim lim E[p T/x (a+(f) a(g))' 
= lim p (a + (f)a(g)) 

L^oo 



(7.3) 



for any T > 0. That is, for any initial state po satisfying the assumptions of the 
theorem, the Boltzmann limit Ft is stationary. This proves Theorem |3.2| □ 



8. Proof of Theorem 13.4 



Theorem 



3.1 



Theorem 3.4 also follows as an almost immediate consequence of our proof of 
Since we are assuming that F G i°°(T 3 ) does not depend on the 



time variable, we have 

U Sl , S2 (u) = exp(z(( S2 - Sl )(£(u) - \(v*F)(u))^ 



(8.1) 



HOMOGENOUS FERMI GAS IN A RANDOM MEDIUM 



o3 



which yields a time-independent shift of the kinetic energy, 

E(p) -> E x (p) := E(p) - X(v*F)(p) . (8.2) 

Thus, under the assumptions that ( |3.17 ) and the crossing estimate (3.18) hold, the 
derivation of the Boltzmann limit reduces to the case treated in [TU] for A = 0, with 
E\(jp) replacing the kinetic energy function E(p). 



In the kinetic scaling limit determined by t 
obtains 



\img[F; m X;T/rf;f,g} 

/j — s-0 



dpf(p)g(p)F T (p) 

for the given choice of F, where Ft satisfies the equation 

d T F T (p) = 2tt / duS(E x (u) - E x {p)) (F T (u) - F T {p)) 



2 and i) — > 0, one accordingly 

(8.3) 



(8.4) 



with initial condition Fo(j>) = F(p). While (8.4| has the form of a linear Boltzmann 
equation, Ft is an auxiliar quantity of which only the stationary solutions are 
relevant for our discus sion . Here, we point out that the renorm alized energy E x in 
the collision kernel in (8.4 1 is determined by F, not by Ft; see (8.2 1. 



Stationary solutions of (8.4 1 are determined by the condition 8tFt(p) = 0, so 
that F T {p) — F (p) = F(p). This holds if and only if F(p) satisfies the self- 
consistency condition 

' f du5(E x {u)-E x (p))F(u), (8.5) 



F(p) 



m x (p) 



where 

m x (p) := 2ir / <h, >){ Exdn-E^p) ) 
This proves Theorem |3.4| 



(8.6) 
□ 
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